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Abstract 

In this paper, a mathematical model of a one-stage Oxford split-Stirling cryocooler has been made to simulate its dynamic 
performance. A set of governing equations has been solved by a numerical method. The present paper emphasizes the validation of 
this model against experimental results for the Oxford split-Stirling cryocooler. The comparison shows acceptable agreement be¬ 
tween the simulation and experimental results. The simulated results on compression work have relative discrepancies of 6.8% (cold- 
tip temperature 80 K) and 6.5% (cold-tip temperature 54 K) with their experimental results respectively. The greatest deviations for 
the minimum pressure and net cooling power are about 6.0% and 8.0% respectively. 

© 2002 Elsevier Science Ltd. All rights reserved. 
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1. Introduction 

Oxford split-Stirling cryocoolers are having been de¬ 
veloped since its invention by the researchers in Oxford 
University at the end of 1970s fl], and many special 
components or technologies such as a linear compres¬ 
sor, no-contact seals, flexure springs, controlling con¬ 
taminants technologies and etc. have been introduced 
during the manufacturing and assembling processes. At 
present, because of their maturity and high reliability, 
Oxford split-Stirling cryocoolers have become quite pop¬ 
ular in many fields, space applications in particular 
[4,5]. 

In an Oxford Stirling cryocooler, the regenerator al¬ 
ways contains a long stack of phosphor bronze gauze 
discs. The compression piston and the displacer move at 
a high frequency (30-60 Hz), and the working gas’s os¬ 
cillating movement is very complicated [3,6,9]. All these 
make it very difficult for the researchers to grasp the 
genuine performance characteristics, and rely on a lot of 
experiments or experimental data during the thermal 
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design process. In order to reduce the development pe¬ 
riod of an Oxford split-Stirling cryocooler, a dynamic 
simulation program has been developed and a visualized 
experimental project has been put forward by us. 

Since Finkelstein’s pioneering the nodal network 
analysis of Stirling Engine [2], a number of nodal anal¬ 
ysis systems are available for the small Stirling cryo¬ 
coolers in the literature [7-9], however, none of them has 
been validated extensively by experiments. To the best 
of the author’s knowledge, a mathematical model for a 
one-stage Oxford split-Stirling cryocooler has been de¬ 
scribed in this paper. A set of differential equations 
based on the conservation of mass, energy, and mo¬ 
mentum for working gas or matrix are solved by a nu¬ 
merical method. The present paper emphasizes the 
validation of the above model against experimental re¬ 
sults for an Oxford split-Stirling cryocooler and ac¬ 
ceptable agreement was fount between the two. 


2. Physical model and governing equations 

The physical model for an Oxford split-Stirling cryo¬ 
cooler is shown in Fig. 1. In order to set up a compara¬ 
tively simplified and practical mathematic model, the 
basic assumptions are made as follows: 
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Nomenclature 


A heat transfer area 

A h heat transfer area for hot chamber 

c p gas specific heat 

c m matrix specific heat 

D h equivalent diameter 

/ friction factor or frequency 

fi, f 2 friction factor 

/ electric current 

m gas mass 

m mass flow rate 

p pressure 

Pr Prantle’s number 

Q heat or cooling power 

R gas constant or electric resistance of winding 
coil 

Re Reynold’s number 

S stroke 

T temperature or period 

u gas flow velocity 

V volume 

W work rate 

x coordinate of flow direction 

Greek letters 

a heat transfer coefficient 

g efficiency or emissivity 

0 differential phase angle 

1 thermal conductivity 

/< dynamic viscidity or modified empirical co¬ 

efficient 
p density 


G 

convergent constant 

T 

time 

y 

ratio of specific heat 

* 

modified resistance coefficient 

Subscripts 

A 

ambient 

C 

compression chamber or compressor 

conduction conduction loss 

cp 

compressor 

D 

displacer 

E 

cold chamber 

effective 

: effective value 

emit 

radiation loss 

f 

boundary of control volume 

gross 

gross cooling power 

h 

hydraulic 

H 

hot chamber or hot end 

i 

the /'th control volume 

input 

compressor input power 

j 

the y'th control volume 

m 

matrix 

max 

maximum 

net 

net cooling power 

out 

outlet 

pump 

pump loss 

R 

regenerator 

reg 

regenerator or regenerative loss 

shuttle 

shuttle loss 

t 

tube or hose 

W 

wall 



Fig. 1. Physical model: (1) compression chamber, (2) connecting hose, 
(3) hot chamber, (4) regenerator, (5) displacer, (6) cold chamber. 

1. one-dimensional (ID) flow and heat transfer; 

2. ideal gas and constant properties; 

3. negligible axial conduction and linear temperature 
distribution of the regenerator matrix; 

4. negligible clearance leakage between the cylinder and 
piston, or the cylinder and displacer; 


5. constant wall temperature of the hot and cold cham¬ 
bers; 

6. momentum and other equations are not coupled. 

The following describe dynamic model and governing 
equations for the five components: namely, compression 
chamber, connecting hose, hot chamber, regenerator, 
and cold chamber, respectively. 

2.1. Adiabatic compressor model and governing equations 

For the one-stage Oxford split-Stirling cryocooler, 
the pressure wave is provided with a linear compressor 
coupled to the hot chamber through a connecting hose, 
as shown in Fig. 1. This design can make the cold-tip not 
be affected greatly by vibration from the compressor. It 
is particularly suitable to coupling with the infrared 
focal plane in many applications. 

Since the compression piston motion is very fast 
and the mixing in the compression chamber is severe, the 
gas temperature and pressure inside the compression 
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chamber are assumed to uniform. The compressor cyl¬ 
inder is assumed to be adiabatic in the analysis, since 
each of the compression and expansion processes occurs 
in such a short period of time that little heat exchange 
between the gas and the cylinder wall could be effectu¬ 
ated. The gas adiabatically compressed in the cylinder is 
assumed to be cooled to room temperature by the cyl¬ 
inder lid and compressor base frame, which were made 
of copper and industry titanium respectively. The cyl¬ 
inder lid and compressor base frame is regarded as the 
after-cooler and assumed to be perfectly effective, so 
that the temperature of the gas leaving it is always equal 
to its wall temperature. 

Applying the first law of thermodynamics to the 
control volume drawn around the volume swept by the 
compression piston in the cylinder, the following equa¬ 
tions can be obtained for ideal gas: 

dfc = — (RT-Amc - — dp c ) (d m c >0) (1) 

Pc V 7 ) 

dVc = V c ( — -y^] (dm c < 0) (2) 

V rn c Pc J 

where m c , V c and p c are, respectively, the mass, volume, 
and pressure of the gas in the control volume; F a is the 
temperature of copper cylinder lid and titanium com¬ 
pressor base frame, which is assumed to be same as 
room temperature; and y is the ratio of specific heats. 
The increment of mass in the above two equations, d m c , 
is related to the mass flow rate into or out of the cyl¬ 
inder, me. If me is defined as positive when the gas is 
leaving the cylinder, the conservation of mass gives 

dm c = ~m c dr (3) 

For the Oxford split-Stirling cryocooler, the com¬ 
pression piston is driven by a linear motor. The equation 
of volume variation for the compression chamber can be 
written as follow: 

Vc=Vc + AV c sin(27t/r) (4) 

where V c and AV C are, respectively, the average volume, 
and volume variation amplitude of the gas in the control 
volume; and / is the frequency for the moving com¬ 
pression piston. In practice, volume variations for the 
compression chamber and cold chamber are not exactly 
sinusoidal wave; so their equations should be rectified 
according to their measured results. 

State equation for the gas in the control volume: 

p c V c = m c RT c (5) 

From the above five equations, the variation with time 
of p c and m c , may be determined by numerical calcu¬ 
lation. The compression power and input power of the 
compressor can be written as follows: 

K = -f <£pcdV c ( 6 ) 


ffinpu, = W cp + J 2 R 


(V) 


where / and R are the current and resistance for winding 
coil of the linear motor. 

The compressor efficiency may be calculated by 


K 

fHnput 


( 8 ) 


2.2. Connecting hose model and governing equations 


The connecting hose links the compressor and the 
expander. The gas flow field is assumed to be ID. For 
governing equations being simplified easily, the second 
order viscous and inertia terms in the momentum 
equation are neglected. However, the second order vis¬ 
cous friction is taking account by a modified resistance 
coefficient t// which is calculated by piecewise linear ap¬ 
proximation [10]. Therefore, the viscous resistance is 
proportional to the mass flow rate with a proportional 
constant i/z depending on the amplitude of the mass flow 
rate. 

Since the connecting hose is small in diameter and 
thick in wall thickness in order to withstand the high gas 
pressure, the gas in the connecting hose is assumed to 
undergo an isothermal process with room temperature 
7],. If we divide the connecting hose into / segments, the 
governing mass and momentum equations for the gas 
in the z'th (i = 1,2,...,/) segment connecting hose are 
given by 


1 8 Qt), 
RTz 0T 

J_ 6(m t ),. 
A t dr 


J_9 

A, 

8 (a),- 

dx 


(m t ) 


0X 

(9) 

o 

II 

J) 

-s-l^ 

+ 

(10) 


As for the modified resistance coefficient t//, the follow¬ 
ing empirical formula for oscillating gas flow in long 
tube from Roach and Bell [11] was adopted. 

i/z = 0.1556[p t (w t ) max d th /7i t ]~°' 201 [(« t ) max M h ] (11) 

where (zz t ) max and d th are the maximum velocity of the 
oscillating gas flows, and the hydraulic diameter of the 
connecting hose respectively. It is worth noting that i jj is 
considered to be a constant during the modeling, but 
should be adjusted by numerical iteration during the 
computation to give a correct value for the responding 
maximum gas velocity (w t ) max . 


2.3. Hot chamber model and governing equations 

Being similar to the compression chamber, the gas 
temperature and pressure are assumed uniform in the 
hot chamber since the mixing effect is severe. The mass 
and energy equations are given by 
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-^=(m t ) 7 -(m r ) in (12) 

OT H c p (Ar H ) AT = oc(^h) t ((Th) t — Tfi w ) (13) 

The volume variation for the hot chamber is in op¬ 
posite phase with that for the cold chamber which was 
given by equation (21). 


compression piston motion. The equations of volume 
variation for the cold chamber can be written as follow: 

V E = V E + AV e sin(2rc/T + <2>) (21) 

Being similar to Eq. (4), Eq. (21) should be rectified 
according to experimental results. 

The gross cooling power of the cooler can be written 
as follow: 


2.4. Regenerator model and governing equations 


The regenerator is divided into many subsystems and 
each subsystem can exchange work, heat and mass with 
the surroundings through its boundaries. The parame¬ 
ters of each subsystem vary with the time and are solved 
with the lumped method. It is defined mass flow rate 
being positive if the gas flows from the hot chamber to 
the cold one and negative for the opposite direction. 
Since the temperature gradient is large within the re¬ 
generator, we divided the regenerator into many sub¬ 
systems. For each subsystem, j (j = 1,2,..., J), the 
governing equations can be written as follows.Conti- 
nuity for gas flow: 

dm . . ,,.. 

— = mj-m j+ i (14) 


Conservation of momentum for the gas (by simplifica¬ 
tion): 


d Pj f 1 Pjtf uj 

dx D h 2 \uj\ 


(15) 


Conservation of energy for the gas: 

C P + Cp(*7Vl “ c V^” lT )j - <Xj A A T Wj 


-^ = 0 


7 ) 


(16) 


Conservation of energy for the matrix: 

0 J 7 

Pm Cm + 0-jAj ( T Wj - Tj) = 0 (17) 

State equation for the gas: 

Pj = p J RT J (18) 


2.5. Cold chamber model and governing equations 


Sgross — / ^ Pe dV E (22) 

We define (^conduction; <2reg; ^shuttle? Qpump* and Q e mit as 
solid conduction loss, regenerator inefficiency loss, 
shuttle loss, pump loss and radiation loss, respectively. 
Other losses of cooling, such as residual gas conduction 
loss, leakage loss, etc. are considered to be negligible. 
Therefore, the net cooling power of the cooler is 


Qnet — Qgross Qe i 


Qrcg ^shuttle Qpnmp Qe i 


(23) 


3. Losses of cooling 

In order to know the net cooling power available at 
the cold end of the cold finger, we need to calculate each 
loss term in Eq. (23). The regenerator inefficiency loss 
and solid conduction loss are the most important. 


3.1. Regenerator inefficiency loss 

In the regenerator of an Oxford Stirling cryocooler, 
the mass flow rate, temperature, and pressure are all 
functions of time. Oscillating flow occurs rather than 
steady state flow. There being considerable void volume 
in the regenerator, the fluctuating gas pressure gave rise 
to the variable mass in the regenerator, and differential 
mass flow rate at different section of the regenerator. In 
view of the fluctuating gas temperature and variable 
mass flow rate, enthalpy flow from hot end to cold end 
of the regenerator, that is to say, regenerative ineffi¬ 
ciency loss can exist. This loss term can be calculated 
through numerical method by 


Being similar to that in the hot chamber, a lumped 
model was applied to the gas in the cold chamber. The 
mass and energy equations are given by 

^=(% )j (I 9 ) 

»zec p (A7e) Ai = 0{(T e ) t ((7e) t — Tew) (20) 

The displacer is driven by a linear motor, and the 
motion phase difference maintains constant ( 1> with the 


2reg =fj> ("i R ) out CpT(Z R , t) dr (24) 

3.2. Solid conduction loss 

The solid conduction loss comprises the heat con¬ 
ductions through the expander cylinder, and the regen¬ 
erator, including the wall (displacer) and the matrix, 
which can be calculated by 
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( ^conduction )r — /Ineffective iP-D^/ 4r) 


Trh — Tv 


(25) 


where /.effective is the effective thermal conductivity to the 
regenerator; and p is a modified empirical coefficient, 
which takes the matrix material, matrix thermal con¬ 
ductivity, matrix stacking, etc. into account. 

Owing to large temperature gradient on the regener¬ 
ator and displacer, the matrix thermal conductivity may 
be varied to a great extent. Therefore, the effective 
thermal conductivity to the regenerator is given by 

l pTrh 

-^effective = ™~— / X{T) dT (26) 

/RH -Ire J Tre 

The heat conduction loss caused by the cold-finger 
cylinder wall can be calculated by 

A f ^ HW 

(2conduction) cy ii n( i er = ~ J 1(T) dT (27) 


can give rise to cooling loss termed as pump loss. The 
pump loss can be calculated by the following equation 
[13]: 


£?shuttle 4.04L d (7h - T e )<5 2 ' 6 

Vc ° CAP). 


nDr 


0.6 


1 1.6 


{Th + T e ) v ^P eak -P eak J 


(33) 


3.5. Radiation loss 

There being large temperature difference between the 
cold-tip and the surroundings with constant ambient 
temperature, radiation loss of cooling power may occur. 
When the coolers works at stable state, the cold-tip 
temperature 7p W is almost constant with a very small 
fluctuation in practice. The radiation loss can be calcu¬ 
lated by 


So the total solid conduction loss is 


Qa 


= (Q 


conduction / r 


(a 


conduction) cylinder 


(28) 


Senut = 5.77 X l(T 4 e ff E 




(34) 


3.3. Shuttle loss 


where 


The displacer/regenerator shuttle heat loss is evalu¬ 
ated from the relation of Baik and Chang [12]: 

^shuttle = uD ° Sd 2n/A D p D (cp) D 7io sin 0 (29) 

where 


Tio — A7 dSd 



(30) 


0 = tan 1 


(ui as) j s/3. T 01 U 2 
(fli + U 2 )/ V2 


(31) 




2d fanf 

h y (Xd 



(32) 


where D d is the displacer outside diameter; S D is the 
displacer stroke; 2 D is the displacer thermal conductiv¬ 
ity; 2 W is the cylinder wall thermal conductivity; p n is 
the displacer density; (C P ) D is the displacer heat capac¬ 
ity; A 7 ’ d is the axial temperature gradient of the dis¬ 
placer; a D and a w are the thermal diffusivity of the 
displacer and the cylinder wall; h is the heat transfer 
coefficient. 


e = 


£ew 


1 


^EW 

^4hw 



(35) 


where e EW and e H w are the emissivities for the cold and 
hot surface respectively; ff EW and ff HW are the area for 
the cold and hot surface respectively. 


4. Friction factor and heat transfer coefficients 

Since we could not find the friction factor and heat 
transfer coefficient for the oscillating gas flow, and the 
transition from laminar to turbulent flow is not obvious 
in the Stirling or pulse tube cryocoolers, we adopted the 
steady friction factor and heat transfer coefficient from 
Ref. [14]. The friction factor is 


/ = max(/i,/ 2 ) (36) 

where the function max is defined as the largest among 
the variables, that is to say, if f\ is equal or larger than 
f 2 , f = fu on the contrary, / = / 2 . f\ and f 2 can be 
calculated by the following formulae respectively 

fx = 1.6+ 180/tfe (37) 


3.4. Pump loss 

Owing to the clearance between the displacer and 
the cold-finger cylinder, some gas is moving periodi¬ 
cally between the clearance and the cold chamber, which 


f 2 = 0.257Re-° 157 ( 38) 

The heat transfer coefficient is calculated by 

a = max(ai,a 2 ) (39) 
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where 


oq = 0.023ARe°' 8 P>'°' 4 /Dh 

(40) 

a 2 = 1.471 ARe°- 5i6 Pi- om /D h 

(41) 

5. Numerical method 



The above set of differential equations governing gas 
flow in connecting hose or regenerator is non-linear and 
unsteady. Some efforts have been made to increase both 
the calculation accuracy and speed. A uniform network 
for space discretization and an implicit scheme for time 
are utilized during the solving process. The control 
volume boundary temperature 7], is given by the second- 


order upwind difference scheme: 


T ti = 1.57%, - 0.57A 2 (m> 0) 

(42) 

7}; = 1.57)+, - 0.57) +2 (m < 0) 

(43) 


The discretization energy equation can be solved by 
the TDMA (tridiagonal matrix algorithm) method. The 
iteration results should be satisfied by periodic stable 
conditions for all parameters and the energy balance 
condition, i.e. the net heat transfer for the matrix of 
regenerator tends to zero during a cycle. The convergent 
equation can be written: 

<f (S0) m = =► 0 (44) 

J t—T 

a can be adjusted according the solution accuracy and 
computational speed. 

6. Computational results and verification by experiments 

Verification of our numerical modeling was carried 
out by comparing the simulation results to experimental 
data. The tested machine is a one-stage Oxford split- 
Stirling cryocooler constructed at Cryogenic Laboratory 
in Shanghai Institute of Technical Physics (SITP). The 
main parameters are as follows: compressor displace¬ 
ment 1.6 cm 3 ; compression piston stroke 8 mm; dis¬ 
placer stroke 4 mm; frequency 40 Hz; initial charge 
pressure 1.3 MPa. The ambient temperature is 303 K. 
Dynamic pressures in the compression and hot cham¬ 
bers were measured with two piezoelectric transducers 
and charge amplifiers. The displacements of the com¬ 
pression piston and the displacer were measured with 
two linear variable-differential position transformers, 
enabling the volume change to be determined accord¬ 
ing to the driven mechanism respectively. The cold-tip 
temperature was measured with a platinum resistance 
thermometer. 

Fig. 2 shows the volume variations versus time for the 
compression and cold chambers. Their measured results 



Fig. 2. Variation of volume with time, <f> = 78°. 


indicate that time for the compression process is as 
about 1.2 times long as that for expansion process in 
some sort; these results deviate the assumption of the 
sinusoidal volume variations for the compression and 
cold chambers. So we utilized equations for the fitted 
curve to experimental results in place of Eqs. (4) and 
(21), respectively. According our experimental results 
while cold-tip temperature maintaining 80 K, the recti¬ 
fied displacement equations can be written by 

V c + AV C sin( 1 Mnf t) (x c ^ 0) ^ 

V c + AF c sin(2.167r/T) (x c < 0) 1 ’ 

= f F E + AK e sin(1 .84 ic/t + 0.92<Z>) (x D ^ 0) 

E \ Ve + Afjj sin(2.167t/T + 1.08$) (xd < 0) 

(46) 

Fig. 3 shows pressure variations versus time for the 
three chambers when the cold-tip temperature maintains 
at 80 K (Fig. 3A) and 54.3 K (Fig. 3B) respectively. The 
simulated dynamic pressures for gas in compression 
chamber, which were derived from Eqs. (4), (21) and 
Eqs. (45), (46) respectively, show that simulated result 
by utilizing rectified displacement Eqs. (45) and (46) has 
higher accuracy to the corresponding experimental re¬ 
sult. In Fig. 3A, the tested pressure amplitudes for the 
compression and hot chambers decrease by 10.8% and 
18.0% of their simulating results respectively; the tested 
pressure amplitude for hot chamber decrease by 12.0% 
of the tested result for the compression chamber. Com¬ 
paring with the simulating result for the compression 
chamber, the simulating results for the hot and cold 
chambers lag 0.3 ms (4.3°) and 0.78 ms (11.2°) in phase 
respectively. The tested results are in phase with the 
simulating results for the compression and hot cham¬ 
bers. The tested maximum pressures for the compression 
and hot chambers are in consistent with the simulating 
ones, however, the measured minimum pressures in¬ 
creasing by about 4.8% and 6.0% than their simulating 
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Fig. 3. Variation of pressure with time, 0 = 78°. (A) Cold-tip temperature 80 K, (B) cold-tip temperature 54 K. 


results respectively. All these pressure amplitudes atte¬ 
nuation and phase delay occur analogously in Fig. 3B 
which the cold-tip temperature maintaining 54 K, how¬ 
ever, the corresponding values are reduced because of 
the increased compressor pressure fluctuation (the sim¬ 
ulated (measured) pressure amplitudes are 307 (283) and 
277 (247) kPa for cold-tip temperature 54 and 80 K, 
respectively), and pressure ratio (the simulated (mea¬ 
sured) pressure ratios are 1.60 (1.54) and 1.46 (1.38) for 
cold-tip temperature 54 and 80 K, respectively). Besides 
all above, Fig. 3A and B indicate that time of com¬ 
pressing process for cold-tip temperature 54 K is longer 
than that for cold-tip temperature 80 K in a working 
cycle. This phenomenon was caused by more load force 
and compressor piston providing with more work on gas 
for cold-tip temperature 54 K than 80 K. By comparing 
Figs. 2 and 3, we can deduce that with the cold-tip 
temperature lowering down, the phase differential angle 
between the volume and pressure variations for gas in 
the cold chamber increases from about 10.2 ms (147°) at 
80 K of cold-tip temperature to 11.5 ms (165°) or so at 
54 K of cold-tip temperature; this increasing phase dif¬ 
ferential angle was the main cause which gave rise to the 
net cooling power rapidly lowering from 630 mW at 80 
K to no net cooling power at 54 K of cold-tip temper¬ 
ature. The above phase lag and amplitude reduction 
may caused by the damp effect of the connecting hose, 
our mathematical model not covering all information 
about the machine structure and movements, and the 
increased dead volume by instillation of transducers. 

Figs. 4 and 5 show the simulated dynamic charac¬ 
teristics of the one-stage Oxford split-Stirling cryocooler 
as a function of run time, including instantaneous 
pressures in three chambers, mass flow rate at the 
compressor outlet, and the hot and cold ends of regen¬ 
erator, and gas temperatures in the compression and 
cold chambers. Fig. 4 shows the mass flow rates at the 
hot and cold ends of regenerator are basically in phase, 



Fig. 4. Variation of mass flow rate and dynamic pressure versus time 
0 = 78°, cold-tip temperature 80 K. 1, 5: outlet of compressor; 2, 4: 
hot end of regenerator; 3, 6: cold end of regenerator. 



Fig. 5. Gas temperature at compression and cold chambers 0 = 78°, 
cold-tip temperature 80 K. 


Temperature (K) 
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and are about 5.57 ms (80.2°) ahead of that at com¬ 
pressor outlet. This is in good agreement with the phase 
differential angle 80° between the compression piston 
and displacer movements. Fig. 5 is a plot of the gas 
temperatures at the compression and cold chambers at 
80 K of cold-tip temperature. These two temperature 
waves are similarly in phase with their pressure waves 
respectively. The compression gas temperature oscillates 
between 305 and 355 K. For 630 mW cooling load at 80 
K of the cold-tip temperature, the cold gas temperature 
fluctuates between 74 and 86 K. 

The dynamic characteristics, including energy rela¬ 
tions, instantaneous pressure, volume variations in the 
three chambers, and P-V diagrams, give a good un¬ 
derstanding of the internal process and performance 
about an Oxford split-Stirling cryocooler. Pressure ver¬ 
sus volume, or P-V plots are very useful for designing 
Stirling cryocoolers. 

Fig. 6 shows comparison of the converged compres¬ 
sor P-V curve with the experimental data. The maxi¬ 
mum pressure predicted agrees well with the measured 
results; however, the minimum pressure is about 4.8% 
(cold-tip temperature 80 K) and 2.3% (cold-tip temper¬ 
ature 54 K) lower. The shapes of both curves are similar, 
and the areas inside the predicted curves are about 17.2 
W (cold-tip temperature 80 K) and 19.7 W (cold-tip 
temperature 54 K) respectively, comparing with that 
inside the measured curves being about 16.1 W (cold-tip 
temperature 80 K) and 18.5 W (cold-tip temperature 54 
K), with relative discrepancies of 6.8% (cold-tip tem¬ 
perature 80 K) and 6.5% respectively. The input power 
of this compressor is approximately 22.1 W with a 
compressor efficiency ( W inpul - I 2 R)/W mput of 78% while 
cold-tip temperature maintaining 80 K; and the input 
power is about 26.8 W with efficiency of 74% while cold- 
tip temperature maintaining 54 K. For cold-tip rated 



Fig. 6. Comparison of simulated P-V diagram against experimental 
results, 0 = 78°. 
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Fig. 7. Cooler performance versus differential phase angle between the 
displacer and compression piston. 


working temperature of 80 K, the total consuming 
power of this tested unit is nearly 38.6 W, including 
approximate 4.5 W for the displacer driven motor and 
about 12 W for the controlling devices. 

Fig. 7 shows the cooler performance versus differen¬ 
tial angle between the displacer and compression piston. 
The compressor input power increases with the increase 
of the phase angle; however, there is an optimized phase 
angle of 78° for the cooling efficiency or 81° for the net 
cooling power respectively. We select operating phase 
angle of 80° in order to harmonize the cooling efficiency 
and net cooling power. 

We simulated the operating cooling performances 
at different cold-tip temperature respectively, as were 
shown in Fig. 8. Fig. 8 indicates that the net cooling 
power is in approximately linear relation with the cold- 
tip temperature. The simulating results are congruent 
with the experimental ones when the cold end temper¬ 
ature is lower than that of liquid nitrogen (77 K); the 
simulating results increasingly deviate the experimental 
ones with increasing cold end temperature in the range 



Fig. 8. Net cooling versus cold-tip temperature, 0 — 78°. 
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Colddown time (min) 

Fig. 9. Cold-tip temperature versus cooling time, <f> = 78°. 

of 80-100 K. When the cold end temperature is 100 K, 
the simulating result is approximately 8.0% larger than 
the experimental one. Fig. 9 shows the cooling process 
for the Oxford split-Stirling cryocooler. This tested 
machine needs about 4.5 min for attaining 80 K with 630 
mW cooling power, and about 10 min for attaining limit 
cooling temperature of approximately 54 K with no 
cooling load respectively. 


7. Conclusion 

The dynamic modeling for a one-stage Oxford split- 
Stirling cryocooler that has been developed can predict 
dynamic process and performances for which the struc¬ 
ture dimensions and operating condition are given. 
Acceptable agreement between experiment and simula¬ 
tion has been found. The developed computer program 
can be utilized to develop an understanding of the in¬ 


ternal process, and predicting the operating perfor¬ 
mances and dynamic characteristics. 
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